N81 - 1471 1 


A THREE-D I MENS I ONAL BODY-FITTED COORDINATE 
SYSTEM FOR FLOW FIELD CALCULATIONS ON 
ASYMMETRIC NOSETIPS 


Darryl W, Hall 
Science Applications, Inc, 
McLean, Virginia 

ABSTRACT 


This presentation describes a three-dimensional body-fitted coordi- 
nate system developed for use in the calculation of inviscid flows over 
ablated, asymmetric reentry vehicle nosetips. Because of the potential 
geometric asymmetries, no standard coordinate system (e.g., spherical, 
axisymmetric reference surface-normal) is capable of being closely aligned 
with the nosetip surface. To generate a 3-D, body-fitted coordinate 
system an analytic mapping procedure is applied that is conformal within 
each meridional plane of the nosetip; these transformations are then 
coupled circumferentially to yield a three-dimensional coordinate system. 
The mappings used are defined in terms of "hinge points", which are 
points selected to approximate the body contours in each meridional 
plane. The selection of appropriate hinge points has been automated to 
facilitate the use of the resulting nosetip flow field code. 
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PROBLEM DEFINITION 


The goal of this effort is the development of a procedure for cal- 
culating supersonic/hypersonic inviscid flows over asymmetric ablated 
reentry vehicle nosetips. These asymmetric shapes, such as illustrated 
in this figure, result from asymmetric transition on the nosetip, which 
occurs at the lower altitudes during reentry (i.e., below 15.24 km). 
Because these shapes occur in the high Reynolds number, turbulent regime, 
with thin boundary layers, an inviscid solution is capable of accurately 
predicting the pressure forces on the nosetip. The nosetip flow field 
solution is also required to provide the required initial data for after- 
body calculations; this coupling of nosetip and afterbody codes allows 
accurate prediction of the effects of the nosetip shape on the afterbody 
flow field. 

The flow field code developed is a finite-difference solution of 
the unsteady Euler equations in "non-conservation" form (i.e., the de- 
pendent variables are the logarithm of pressure, P, the velocity compo- 
nents, u,v,w, and the entropy, s). In this approach the steady flow 
solution is sought as the asymptotic limit of an unsteady flow, starting 
from an assumed initial flow field. 


CALCULATION OF SUPERSONIC/HYPERSONIC INVISCID FLOWS OVER ASYMMETRIC ABLATED 
REENTRY VEHICLE NOSETIPS 



APPROACH 

• FINITE-DIFFERENCE SOLUTION OF UNSTEADY EULER EQUATIONS 

• STEADY FLOW SOLUTION SOUGHT AS THE ASYMPTOTIC LIMIT OF 


UNSTEADY FLOW 


COORDINATE SYSTEM REQUIREMENTS 


It is well known that accurate numerical calculation of fluid flows 
requires the use of a coordinate system closely aligned with the principal 
features of the flow. For the nosetip problem this requirement would be 
satisfied by a coordinate system which closely follows the body shape and, 
hence, the streamlines of the flow. Because of the asymmetric nosetip 
geometries being considered, standard coordinate systems (e.g., spherical, 
axisymmetric reference surface-normal) are incapable of being aligned with 
the nosetip surface at all points. Thus, a coordinate transformation is 
sought that will align the coordinate system with an arbitrary nosetip 
geometry. By requiring the transformation to be in analytic form, the 
need of solving partial differential equations to define the transformation 
can be avoided. Finally, the transformation should be in a form that 
readily lends itself to automated definition, minimizing the inputs re- 
quired of a user of the code. 


OPTIMUM COORDINATE SYSTEM FOR NUMERICAL FLOW FIELD CALCULATIONS 
IS BODY-ORIENTED 


COORDINATE TRANSFORMATION SOUGHT THAT: 

1.) ALIGNS COORDINATE SURFACES WITH BODY 
SURFACE 


2.) IS ANALYTIC (SOLUTION OF PDE’S NOT REQUIRED 
TO DEFINE TRANSFORMATION) 


3.) CAN BE READILY AUTOMATED (TO MINIMIZE INPUTS 
REQUIRED FROM USER) 
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COORDINATE TRANSFORMATION 


V. 


The nosetip geometry is defined in an (x,y,4>) cylindrical coordinate 
system, and a mapping to a (£»n,0) transformed coordinate system is sought. 
Since current reentry vehicle nosetips are initially axisynmetric (prior 
to ablative shape change), it is assumed that nosetip cross-sections re- 
tain some "axi symmetric" character during reentry. Thus, no transforma- 
tion of the circumferential coordinate is required, and 0 = $ is assigned. 
(This transformation can readily be generalized to 0 = f (<p ) if required 
for other applications of this approach.) Within a f = constant merid- 
ional plane, the transformation reduces to the two-dimensional form 
S = £(x,y), n = n(x,y). Conformal transformations from the z = x+iy to 
thec= ?+in plane are desirable, ensuring that an orthogonal (£,n) grid 
maps back onto an orthogonal grid in the (x,y) plane. 


(x,y,<f>) CYLINDRICAL COORDINATES IN PHYSICAL SPACE 


(C.n.e) COORDINATES IN TRANSFORMED SPACE 


TRANSFORMATION OF CIRCUMFERENTIAL COORDINATE NOT REQUIRED 
(NOSETIPS INITIALLY AX I SYMMETRIC ) ; ASSUME TRANSFORMATION 
TAKES THE FORM 

5 — E (x ,y ><p ) 

n =n (x,y,<fi) 

e = 


IN A MERIDIONAL PLANE (<p = CONSTANT), THE TRANSFORMATION 
REDUCES TO 


£ = £ (x.y) 
n = n (x.y) 


V. 
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DEFINITION OF TRANSFORMATION 


The approach used to define the coordinate transformations is a 
modification of the "hinge point" approach of Moretti*. The mapping 
is defined as a sequence of conformal transformations of the form 

z j + i - 1 = [z j - h j+i,/ J 

where zj = Xj +iyi (j = 1 is physical space) and hi is the ith hinge 
point in the Zj plane. The hinge points in the physical (zi) plane are 
selected to approximately model the body geometry. By mapping the 
hinge points sequentially onto the horizontal axis, the image of the 
body surface will then be a nearly horizontal contour. 

INDEPENDENTLY IN EACH MERIDIONAL PLANE, DEFINE A SEQUENCE 

OF CONFORMAL TRANSFORMATIONS 


Z j + 1 ■ 1 “ lz j - W J = 


z. = x. + 1 j (j = 1 IS PHYSICAL SPACE) 


h < . = iit "HINGE POINT" IN j— SPACE 

» »J 

HINGE POINTS ARE SELECTED TO APPROXIMATE BODY GEOMETRY 



*Moretti, G. , "Conformal Mappings for Computations of Steady, Three- 
Dimensional, Supersonic Flows," Numerical /Laboratory Computer Methods 
in Fluid Mechanics, ASME, 1976. 


SEQUENCE OF TRANSFORMATIONS 


In the jl!l mapping of the sequence, the transformation is centered 
around the hinge point hj+jj. The mappings have the property of keep- 
ing the hinge points, h-jj (i< j+1) on the horizontal axis, while mapping 
the hinge point hj+ 2 ,j onto the horizontal axis. Thus, after JA trans- 
formations, all ja+ 2* hinge points in the JA+1 space will lie on the 
horizontal axis. (Each mapping in this sequence may be considered a 
"point-wise Schwa rz-Christoffel " transformation.) This figure illustrates 
the sequence of transformations for JA = 3. 



EXPONENTS OF TRANSFORMATIONS: 6 - = 
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TRANSFORMATIONS - CONTINUED 


In order to establish a grid suitable for flow field calculations 
when the image of the body contour is a nearly horizontal surface, it is 
desirable to have the image of the centerline external to the body lie 
along the vertical axis. This is achieved using an additional conformal 
transformation, centered around the second hinge point, of the form 

1/2 

z JA+2 = ^ Z JA+1 ' h 2,JA+l^ 

The last transformation is a simple stretching (which is also conformal): 


C Z + in - azj /\ + 2 . 


(This stretching is used in the calculation procedure along the center- 
line.) This figure illustrates the body contour resulting in the £-plane 
for the case of a sphere with JA = 3, where the body surface is defined 

as n = b(£). 


MAP CENTERLINE ONTO VERTICAL AXIS WITH 

z = (z . h I 1/2 
z JA+2 Z JA+1 n 2,JA+r 

ALLOW FOR SIMPLE STRETCHING (REQUIRED FOR CENTERLINE 
TREATMENT) WITH 

C = 5 + i n = az j A +2 


RESULTING BODY CONTOUR: 
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COMPUTATIONAL TRANSFORMATION 


For the flow field calculation it is desirable to have equally 
spaced grid points. Thus, a transformation to a computational coordinate 
system (X,Y,Z) is used, in which grid points are equally spaced circum- 
ferentially in 6, longitudinally in £ within each meridional plane, and 
in n between the body and the shock. It is important to note that the 
(X,Y,Z) system is not orthogonal, and that the computational transforma- 
tion varies with time as the bow shock position varies during the time- 
dependent calculation. These sketches illustrate the computational grids 
resulting in a meridional plane in both physical (z = x+iy) and trans- 
formed (? = £+in) space for a typical ablated nosetip contour (with the 
shock layer thickness exaggerated for clarity). 


DESIRE GRID POINTS EQUALLY SPACED IN Q ALONG BODY. IN n BETWEEN 
BODY AND SHOCK, AND IN 6 CIRCUMFERENTIALLY 



v = l 7 . n-b(C. e) 

IlToT * cU,e.T)-b(c,ey 
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PARAMETERS OF THE TRANSFORMATION 


In transforming the governing equations from physical to the (X,Y,Z) 
computational coordinates, certain derivatives of the transformation are 
required. Because the transformation has been defined in analytic form, 
these derivatives can readily be evaluated analytically and are functions 
only of the hinge point locations. Within a meridional plane (<f> = constant), 
the required derivatives are g = 3c/3z and <!> = 3(log g)/3£. Circumferen- 
tially, the independent transformations in each meridional plane can be 
coupled to produce a three-dimensional transformation by assuming that 
hinge point locations can be expressed as hi The required circumfer- 
ential parameters of the transformation, z,& and g<j,, can be evaluated 
analytically if each meridional plane has the same number of hinge points 
and assuming the form of interpolating functions for hi,j(<J>). Alterna- 
tively, it has been found to be sufficient to evaluate and g^ from 
Taylor series expansions using data at computational (X,Y,Z) mesh points, 
with the forms of the resulting expressions shown in the figure. 


REQUIRED IN WRITING GOVERNING EQUATIONS IN TRANSFORMED 
COORDINATES 


9 “ H = Ge1 “ ' + S = 'S + 9 y 

* = iO°3_3l 
3C 

CAN BE EVALUATED ANALYTICALLY 

CIRCUMFERENTIAL PARAMETERS OF THE TRANSFORMATION 

CAN BE EVALUATED ANALYTICALLY IF EACH MERIDIONAL 
4* 

PLANE HAS THE SAME NUMBER OF HINGE POINTS, ASSUMING INTER- 
POLATING FUNCTIONS FOR h. .{<>) 

i » J 

ALTERNATIVELY, EVALUATE FROM, TAYLOR SERIES EXPANSIONS: 

C ? -Ci-9(z z -z,) 

C, = E.+ In* = - . J - L - 

<P 4>2 -( Px 

92-9J-9 2 , * l (z 2 -Z 1 ) 

g = 

( ■+ (X-AX.Y.Z), ( ) 2 - (X+AX,Y,Z) IN COMPUTATIONAL MESH 
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AUTOMATIC GENERATION OF HINGE POINTS 


To simplify the application of this coordinate transformation to the 
asymmetric nosetip flow field problem, the selection of hinge points that 
define the transformations has been automated. Within each meridional 
plane to be computed, body normals are constructed at points equally 
spaced in wetted length along the body profile. The hinge points are then 
selected to lie a distance 6 inside the body along these normals. By re- 
lating 6 to any convenient scale factor for a nosetip geometry, the only 
input required of the user of the code is the number of hinge points to 
be used. The locations of the first two hinge points (i.e., those that 
lie on the x axis) are the same in each meridional plane, in order to 
simplify the treatment of the centerline. Typically, no more than nine 
hinge points per meridional plane (JA = 7) are necessary for the nosetip 
flow field problem. 



HINGE POINTS LOCATED DISTANCE 6 ALONG INWARD BODY NORMALS, FROM 
BODY POINTS EQUALLY SPACED IN WETTED LENGTH 


ONLY INPUT REQUIRED OF USER IS NUMBER OF HINGE POINTS TO BE 
USED IN EACH MERIDIONAL PLANE 
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TREATMENT OF CENTERLINE 


The greatest complication encountered in the use of this 3-D coordi- 
nate transformation is the extra care that must be taken in treating the 
grid points on the centerline. Since the transformations in each merid- 
ional plane are independent, the scale factors g = 3;/3z along the 
centerline will not be the same in each meridional plane. Thus, one 
computational grid point at the centerline will represent different physical 
points for each value of <t> . To minimize these discrepancies, the stretch- 
ing transformation x, = azj /\+2 is used to ensure that the images of the 
first hinge point are coincident in all meridional planes. The remaining 
discrepancies are small enough that simple linear interpolations can be 
used to account for differences in the scale factors. 

In addition to the mapping complications along the centerline, the 
governing equations in cylindrical coordinates are singular along y = 0. 

This difficulty has been avoided by using a Cartesian (x^,x 2 ,X 3 ) coordinate 
system for the centerline analysis. The required Cartesian derivatives 
can be expressed in terms of the radial derivative 8/8y in cylindrical 
coordinates for certain values of <j>, as shown in this figure. The only 
restriction resulting from this analysis is that computational planes must 
be located at <t> = 0, m/2, it, and 3tt/2. 

AT THE CENTERLINE (y * 0), SCALE FACTORS (9 = 3 ^/3z } VARY WITH * 


STRETCHING TRANSFORMATION USED TO MINIMIZE DISCREPANCIES, WITH 


1 „ V .JA+2^ = 0) 
k h l,JA+2^ * 


CARTESIAN COORDINATES (xj.x^Xj) USED IN CENTERLINE ANALYSIS 


3x j 9x 


J_ = cos$i-- sin* X 


9Xr> 


ay y 3<J> 


3 - 3 cos 3 3 

3^ - sin *3y + T* 


WITH 


Tim \ 3 x r 

y+0 y ckjT 9y9<f> 


FINITE, 


db ■ cos * h • * = 0> * 


3 _ cin . 3 , . tt 3rr 

3x, ' s,n * 3y ’ * * 2 ’ 2 


325 


RESULTING FLOW FIELD CODE 


The 3-D, time-dependent, inviscid nosetip flow field code that was 
developed using the 3-D coordinate transformation described here is called 
CM3DT (£onformal Mapping 3-jl X ransoni ' c ) • ' This coc ^ e can treat ideal or 
equilibrium real gas thermodynamics, has both pitch and yaw capability, 
and is able to treat weak embedded shocks on indented nosetips using the 
^-differencing scheme*. To provide total body inviscid flow field capa- 
bility, the CM3DT code has been coupled to the BM0/3IS**, NSWC/D3CS$ + , 
and STEIN ++ afterbody codes. Complete details on the CM3DT analysis and 
results obtained with this code may be found in the following references: 

Hall, D. W., "Inviscid Aerodynamic Predictions for Ballistic Reentry 
Vehicles with Ablated Nosetips," Ph.D. Dissertation University of Penn- 
sylvania, 1979. 

Hall, D. W. , "Calculation of Inviscid Supersonic Flow over Ablated Nose- 
tips," AIAA Paper 79-0342, January 1979. 

CHjr.T 'CONFORMAL MAPPING 2-D TRANSONIC} 

N05FTIP f LOW FIELD LODE 


• I PEAL Oft EQUILIBRIA REAL GAS THERMODYNAMICS 

• PITCH AND YAW CAPABILITY 

■ .V DIFFERENCING scheme LTSEO TO treat WE At'. EMBEDDED 
SHOCKS ON INDENTED NOSETIPS 

• COuPLFD TO AFTfRBCOY COOES FOR T0rA| JNVTSCID 
FLOW FIELD CAPABILITY 

• BMO/3IS 

• ISSWC/D3C5S 
« ST? Ih 

• 8I,OOC 10 CORE STORAGE REQUIRED 


*Moretti, G., "An Old Integration Scheme for Compressible Flow Revisited, 
Refurbished, and Put to Work," Polytechnic Institute of New York, POLY- 
M/AE Report 78-22, September 1978. 

**Kyriss, C. L. and Harris, T. B. , "A Three-Dimensional Flow Field Computer 
Program for Maneuvering and Ballistic Reentry Vehicles," 10th U. S. Navy 
Symposium on Aeroball istics, July 1975; also, Daywitt, J., Brant, D., and 
Bosworth, F. , "Computational Technique for Three-Dimensional Inviscid Flow 
Fields about Reentry Vehicles, Volume I: Numerical Analysis," SAMSO TR- 

79-5, April 1978. 

+ Solomon, J. M. , Ciment, M., Ferguson, R. E., Bell, J. B., and Wardlaw, 

A. B., Jr., "A Program for Computing Steady Inviscid Three-Dimensional 
Supersonic Flow on Reentry Vehicles, Volume I: Analysis and Programming," 

Naval Surface Weapons Center, NSWC/WOL/TR 77-28, February 1977. 

++ Marconi, F., Salas, M., and Yaeger, L., "Development of a Computer Code 
for Calculating the Steady Super/Hypersonic Inviscid Flow around Real 
Configurations, Volume I. Computational Technique," NASA CR-2675, April 1976. 
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CM3DT RESULTS 


This figure presents some typical results obtained with the CM3DT 
inviscid nosetip flow field code. Shown are comparisons of predictions 
to data obtained for the PANT Triconic shape* at M =5. It is signifi- 
cant that attempts to compute the flow over this slender shape using a 
time-dependent code formulated in a spherical coordinate system were 
unsuccessful. CM3DT, with its body-oriented coordinate system, was able 
to obtain converged solutions for this shape, with the predictions agree- 
ing well with the data, as seen in this figure. 





*Abbett, M. J. and Davis, J. E., "Interim Report, Passive Nosetip Tech- 
nology (PANT) Program, Volume IV. Heat Transfer and Pressure Distri- 
bution on Ablated Shapes, Part II. Data Correlation and Analysis," 
Space and Missile Systems Organization, TR-74-86, January 1974. 
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CM3DT - RUN TIMES 


On a CDC Cyber 176 computer, the CM3DT invlscid nosetip code with 
A-differencing requires approximately 0.00045 CP seconds per grid point 
per time step (iteration). Typically, 400-500 time steps are required 
to obtain a converged solution. It is estimated that the computer time 
required for a solution has been increased by approximately 20% by using 
the 3-D coordinate transformation described here, when the parameters 
of the transformation on the moving grid are updated every ten time steps. 
When compared to the standard MacCormack differencing scheme, the use of 
A-differencing scheme increases the run time requirements approximately 
50% for this code. 


ON A CDC CYBER 176, CM3DT REQUIRES 0.00045 CP SECS/POINT/STEP 
FOR IDEAL GAS CALCULATIONS WITH A-DIFFERENCING 


t 20% PENALTY INCURRED FOR COORDINATE 
TRANSFORMATION (PARAMETERS ON MOVING 
GRID UPDATED EVERY 10 TIME STEPS) 

t 50% PENALTY INCURRED FOR A-DIFFERENCING 
(RELATIVE TO MAC CORMACK DIFFERENCING) 
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